Lanczos study of the S 



1/2 frustrated square-lattice antiferromagnet in a magnetic 
field 
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We study the zero-temperature phase diagram of the frustrated square-lattice S — 1/2 antiferro- 
magnet in an external magnetic field numerically with the Lanczos algorithm. For strong frustration, 
we find disordered phases at high (and low) magnetic fields. Between these two disordered phases we 
find a plateau in the magnetization curve at half of the saturation magnetization which corresponds 
to a state with up-up-up-down (uuud) spin order. This and other considerations |1| suggest an 
unusual ordering scenario: There are an ordered phase with a spin gap (the plateau) and disordered 
magnetically gapless phases above and below. 

The transition to saturation is studied in further detail and problematic conclusions in earlier in- 
vestigations of this region are pointed out. 



PACS numbers; 75.50.Ee, 75.40.Mg, 75.45. -Fj 
I. INTRODUCTION 

Frustrated magnets are known to exhibit rich physics 
and diverse properties in zero external field (compare this 
volume and B]). Exciting behavior in an external mag- 
netic field was also observed recently in several cases: For 
instance, the S = 7/2 three-dimensional antiferromagnet 
Gd3Ga50i2 is disordered (but gapless) in zero external 
field and orders only in finite magnetic fields |H]. An- 
other example is the recently discovered two-dimensional 
S = 1/2 spin-gap material SrCu2 (603)2 which exhibits 
several plateaux in the magnetization curve B . The ori- 
gin of these plateaux and the nature of the corresponding 
spin states are under intense debate Q and it is suggested 
that some of the plateaux in SrCu2 (603)2 might be de- 
scribed by certain types of ordered states. 

A plateau in the magnetization curve of the triangu- 
lar lattice at one third of the saturation nragnetization 
is theoretically well known (see, e.g., ^). However, gen- 
eral magnetization plateaus in two dimensions (including 
those of SrCu2 (603)2) are still not well understood de- 
spite recent progress on the subject M. 



II. THE MODEL 

In PI we have studied the aforementioned issues in one 
of the simplest two-dimensional frustrated magnets, the 
frustrated square lattice antiferromagnet whose Hamilto- 
nian is given by 



H 
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where the first sum is over nearest neighbor pairs (x, y) 
and the second one over diagonal neighbor pairs [x, y\ . 
Here we will present details of m and further results ob- 
tained from a numerical diagonalization of the S = 1/2 
version of the Hamiltonian (1^) with the Lanczos method. 



At zero field, the frustrated square lattice S = 1/2 
Heisenberg model has been studied extensively using also 
exact diagonalization p-p^. One reason for the popu- 
larity of this model is that it exhibits melting of long- 
range magnetic order by quantum fluctuations: For small 
frustration J' ^ J, one finds antiferromagnetic Neel or- 
der while for J' ^ J a certain type of coUinear order 
described by a wave vector (tt, 0) or (0,7r). Finally, in 
an intermediate region around J' = J/2, a spin-liquid 
groundstate with a gap to magnetic excitations is found. 

In the presence of an external field, an important quan- 
tity is the magnetization 



W-^fe^l' 



SV 
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which we normalize to saturation values | (M) | = 1 . Mag- 
netization curves for the model ( p]) h ave already been 
computed by exact diagonalization ]ll] , p^ . However, the 
study of [p] was restricted to a special lattice of 4 x 4 sites 
which can also be interpreted in terms of several other 
geometries. For instance, a strip of width 4 is equivalent 
to a frustrated four-leg ladder and then plateaux are ex- 
pected in the magnetization curve at (M) = 0, 1/4, 1/2 
and 3/4 (see, e.g., Q). Precisely these values were also 
observed in ||ll|, but it remains to be clarified whether 
they are an artifact of the special geometry. The other 
study ||l^ excluded the region J' > J/2 which we find 
to be the most interesting one and boundary conditions 
were used in |12| which frustrate the order which wants 
to develop. 



III. SINGLE MAGNON DISPERSION 

It is instructive to look first at the dispersion of a single 
magnon above the ferromagnetic state with (M) = 1 
(all spins aligned along the field). When we formally set 



h = Q, the one-magnon dispersion for the S = 1/2 model 
is given by 



£is{kx,ky) = — 2(J + J') + J {coskx + cosky) 
+2 J' cos kx cos ky . 



(3) 



For J' < J/2, this dispersion has a single minimum at 
kx — ky = TT corresponding to Neel order. On the other 
hand, for J' > J/2 two equivalent minima are found at 

which signal coUinear 



U, h^y 



TT OT kx — TT, k,i 



order. The case J' = J/2 is special: Here lines of minima 
are found for either kx = n or ky = n . 



IV. MAGNETIZATION CURVES 
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FIG. 1. Zero-temperature magnetization curves of the 
S = 1/2 model at J' /J = 0.6. The lines are for clusters of 
size 4x4 (dashed), 6x4 (dotted), 6x6 (full) and 8x6 
(dashed-dotted) (Inset: 8x8 (dotted) and 10 x 10 (dashed)). 
The bold line is an extrapolation (compare the text for de- 
tails). 

We now discuss the behavior in the entire field range 
using exact diagonalization. Guided by the dispersion 
(0), we have considered only rectangular clusters with 
an even number of sites in both directions. Finite-size 
magnetization curves are shown in Fig. ^for J'/ J — 0.6. 
As we have pointed out in section ||, finite-size effects 
must be examined carefully. Even if the (M) = 3/4 and 
1/4 plateaux are quite pronounced on the 4x4 and 6x4 
clusters, we believe that they are artifacts of the special 
system sizes. Indeed, the data for wider strips indicates 
a vanishing plateau at (A/) = 3/4 in the thermodynamic 
limit. In contrast, the plateau at (A/) — 1/2 is only a 
little narrower on the 6x6 cluster than for the smaller 
clusters. In the extrapolation (bold line in Fig. nl), we 
have therefore drawn a plateau with (M) = 1/2 but none 
at (M) = 1/4 and 3/4. The value of the spin gap {i.e. the 
boundary of the (Af ) = plateau) has been taken from 
WM. Otherwise the extrapolation was obtained by the 
standard procedure of connecting the mid-points of the 



steps of the magnetization curves at the largest available 
system size. 

Fig. shows the projection of the magnetization curves 
onto the horizontal axis, i.e. the locations of the finite- 
size jumps in the magnetization curves. From this one 
can then first read off how the width of the (A/) = 1/2 
plateau varies with J' / J . An analysis of its finite-size 
behavior |1| indicates that the {M) — 1/2 plateau exists 
in the region 0.5 < J' / J < 0.65. 
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FIG. 2. Transition lines for the S = 1/2 model. The lines 
are for clusters of size 4x4 (dotted), 6x4 (dashed) and 6x6 
(fuU). 

Fig. H further shows that the transition to saturation 
is special at J' — J/2, as is expected from (H): A jump of 
size 5{M) — 1/L occurs just below saturation where L is 
the length of the shorter edge of the cluster. This implies 
a smooth transition to saturation in the thermodynamic 
limit. The particularity of the transition to saturation 
for J' « J/2 was already noted in [|2| where, however, 
the jump was obscured by the choice of cluster geome- 
tries. We therefore believe the asymptotic form proposed 
in ||l^ just to be a manifestation of these boundary ef- 
fects in combination with the questionable assumption 
that h{{M)) has a power series expansion. It should be 
possible to obtain the correct asymptotic form from Bose 
condensation of magnons p5[ into the lines of minima of 
the single magnon dispersion (p|), but this computation 
has not yet been carried out. 

For J' « J/2, the magnetization curves further show a 
pronounced finite-size plateau at (Af ) = 1 — l/L, i.e. just 
below the jump. It was actually proposed in pq | that a 
plateau just below saturation should exist for J' > J/2, 
but our candidate disappears as L — > c» and we sus- 
pect an error in the treatment |T3] of the single-magnon 
dispersion (||). 



V. STATIC STRUCTURE FACTORS 

To gain further insight into the phase diagram and 
in particular in order to understand the nature of the 



(M) = 1/2 plateau state better, we use the static struc- 
ture factors which are defined asQ 






The finite magnetization leads to a trivial peak in 
S'^^(0, 0). We observe additional peaks in S°'^{kx,ky) at 
(0,7r) or (tTjO) and (tt, tt). The evolution of these peaks 
with J' /J is shown in Figs. and H for two values of (M). 



(M) = 1/2 and for 0.5 < J' /J < 0.65, the longitudinal 
components develop uuud order with three spins point- 
ing up and one down on a four-spin plaquette. While it 
may be possible that some order exists in S'^^(0,7r) at 
(M) = 2/3 for intermediate J' (see Fig. ^), the structure 
factors at other magnetizations and Monte Carlo simula- 
tions of the classical model indicate disordered states for 
magnetizations sufficiently far above or below the uuud 
state at (M) = 1/2 |]. 
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FIG. 3. Peaks of the static structure factors of the S — 1/2 
model as a function of J' /J at (M) = 2/3 on the 6x6 clus- 
ter. The longitudinal structure factor 5"'^(0,7r) is shown by 
triangles pointing up, S'''' (77,71) by triangles pointing down 
and the transverse structure factor S"^^(0,7r) by circles and 
S'^^(7r,7r) by squares. 
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FIG. 4. Same as Fig. || but for (Ad) = 1/2 and on the 
6x6 cluster (filled symbols) as well as on the 4x4 cluster 
(open symbols). 

Neel and collinear order develops in the transverse 
components for small and large J' , respectively. Just at 



* The symmetrized variant (S°''^{k^,ky) + S"'^{ky,k^)) /2 
is used when the cluster is not square or the groundstate not 
symmetric under reflection across the diagonal of the cluster 
{x ^ y). 



VI. GROUNDSTATE QUANTUM NUMBERS 

For an odd number of magnons V/2 — S^ , one expects 
the groundstate to carry k^ = fcy = tt (corresponding to 
Neel order) for J' <^ J while for J' ^ J one expects 
collinear order and thus kx = 0, ky — tt or k^ = tt, 
ky = 0- compare the discussion of the one-magnon sector 
in section [IF When V/2 — S^ is even, the groundstate 



carries momentum kx — ky = 0. However, as was pointed 
out in |||,|l2l , groundstate level crossings are still expected 
in the sectors with V/2 — S^ = 2n and n odd. Then the 
groundstate is even under diagonal reflections for J' <^ J 
and odd for J' 3> J. 
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FIG. 5. Positions of unique groundstate level crossings. 

Triangles are for the 4x4 cluster, squares for the 6x6 cluster 

and circles for the 8x8 cluster. 

The unique groundstate level crossings are shown in 
Fig. H (in some sectors more than one groundstate level 
crossing is observed). The crossing at (M) = on the 
6x6 cluster was taken from [nOJ. Note that according 
to (||), the crossing point must move to J' = J/2 for 
(M) -^ 1 where a first-order transition between Neel or- 
der and collinear order occurs. The position of the cross- 
ing at (M) = 1/2 agrees with our earlier estimates for the 
transition between the uuud plateau state and collinear 
order. Also the crossing at {M) = is close to recent 
estimates (see, e.g., |1J]) for the location of the zero-field 
transition between the disordered spin-liquid state and 
the collinear state. 



VII. SUMMARY 

We have studied the frustrated square lattice S — 1/2 
Heisenberg model in the presence of a strong magnetic 
field using the Lanczos method. A plateau with (Af ) ~ 
1/2 exists in the region 0.5 ^ J' /J ^ 0.65 and corre- 
sponds to an up-up-up-down ordered state. The trans- 
verse spin components show Neel and collinear order for 
small and large J', respectively. The intermediate region 
contains the ordered (A/) — 1/2 state and disordered re- 
gions at higher and lower fields ||] . Finally, we have found 
that the transition to saturation is special at J' = J/2. 
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